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Abstract The main features of the conductivity of doped single layer graphene are 
analyzed, and models for different scattering mechanisms are presented. Many possible 
dependencies of the cross section on the Fermi wavelength are identified, depending on 
the type of scattering mechanism. Defects with internal structure, such as ripples, show 
non monotonous dependencies, with maxima when the Fermi wavelength is comparable 
to the typical scale of the defect. 



1 Introduction. 

A great research effort devoted to graphene started after the realization that single 
layer graphene can be isolated and that the number of carriers can be tuned [T][2]. The 
properties of epitaxially grown samples with few graphene layers [3], has also induced 
a significant activity. Many of the potential applications are related to the design of 
electronic devices. 

We analyze here some topics related to the electronic transport properties of single 
layer graphene. We do not consider models for the carrier mobility in systems with more 
than one layer. It is interesting to note in this respect that the mobility of graphite is 
significantly higher than that of few layer graphene samples 2]. We will also not analyze 
the interesting topic of the transport properties of undoped graphene. 

We restrict ourselves here to doped single layer samples in the diffusive regime, 
where the mean free path is shorter than the sample size. In this regime, the motion of 
the carriers can be described semiclassically, and the carrier mobility is determined by 
the different scattering mechanisms present in the system. If we describe the effects of 
the scattering events on the carrier distribution function using the Boltzmann equation, 
which is a reasonable approximation in this limit, the conductivity of graphene is given 
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by: 

a = iV(,p)^(,p) (1) 

where A''(ep) is the density of states, and T{e-p) is the scattering time. 

We present, in the following section, a brief description of the behavior of the con- 
ductivity in single layer graphene, as function of carrier concentration and temperature. 
Then, we analyze different scattering mechanisms which may limit the mobility at low 
temperatures. 

2 Qualitative behavior of the conductivity in graphene. 

2.1 Low temperature limit. 

It was realized at an early stage that the conductivity at low temperatures was almost 
directly proportional to the carrier density [1], and this dependence has been repeat- 
edly confirmed[5]- As, N{e-p) oc in graphene, where p is the carrier density, this 
behavior implies that T{e-p) oc y^p. On the other had, for weak local scatterers the 
Born approximation predicts an inverse scattering time proportional to the density of 
states, oc ni^pN {e-p) oc ^/p, where rij^p is the number of scatterers. Inserting 
this expression in eq.([T}, we find that weak scatterers lead to a conductivity which is 
independent of the carrier density, in contradiction with the observed behavior. 

The linear dependence on carrier density of the conductivity implies that the scat- 
tering time should increase with density as r(p) oc y^, or, analogously, the cross section 
of the defects, which in two dimensions is given by a length, should scale as fcp ^. 

The first mechanism proposed compatible with a t oc ^/p dependence was scat- 
tering by charged impurities [6l[7]- The Coulomb potential is scale invariant, as it only 
depends on the product of the charge of the impurity and the electron charge, Ze^ , 
which can be rendered dimensionless by dividing it by the Fermi velocity. Hence, the 
cross section should be proportional to the only length scale in the problem, which 
is the Fermi wavelength, Ap = fcp^. This dimensional argument [8] remains valid even 
when screening by the carriers is included, as the Fermi-Thomas screening length is 
proportional to the Fermi wavelength in graphene. The electronic structure of graphene 
near a charged impurity has been studied extensively in later times [Ol ITOllllj . without 
changing the previous analysis. 

A different mechanism which leads to a linear dependence on the Fermi wave- 
length of the scattering cross section is induced by strong scatterers, such as lattice 
vacancies [12]. These defects change significantly the local density of states, as they 
induce partially localized states at the Dirac energy. As reviewed below, the scattering 
phasehifts cannot be described by the Born approximation [T3lll4lll5] . 

Other mechanisms which lead to deviations from the T~^{fp) oc N{ep) dependence 
expected from the Born approximation are scattering by ripples [16] (see also below), 
and scattering by defects which lead to long range distortions of the lattice, such as 
dislocations. 

2.2 Finite temperature conductivity 



The resistivity of single layer graphene rises as function of temperature [17|[T8iri9] . and 
deviates significantly from its low temperature value at room temperature. 
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The coupling to in plane phonons is well understood [20ll2T| . The phonon band 
width is ~ 0.2eV. The number of thermally excited phonons at room temperature is 
too small to explain the observed rise in the resistivity. Single layer graphene can also 
support out of plane fiexural modes, which show a high density of states at low energies. 
The coupling to the electrons is quadratic on the phonon coordinates, and the resulting 
scattering rate is also too low to explain the observed temperature dependence of the 
resistivity [22]. Note that these modes can be pinned by the substrate [23], reducing 
their density of states at low energies. 

The most likely explanation for the temperature dependence of the resistivity is 
scattering by substrate modes i 24 j. The most common substrate used in experiments 
on graphene samples obtained by mechanical cleavage is Si02, which is a polar insu- 
lator. The electrons in graphene couple to the electric fields induced by the surface 
polar modes. These modes can be thermally excited at room temperature. A fit using 
the observed frequencies and dielectric function of Si02 gives a good agreement with 
experimental data[19j. 



2.3 Effects of the substrate 

Localized charges in the substrate lead to Coulomb scattering and modify the low tem- 
perature mobility. At finite temperatures, the conductivity depends on the coupling of 
the carriers to the substrate modes. Hence, the transport properties of single layer 
graphene are determined, to a large extent, by the substrate and by the general prop- 
erties of the surrounding environment. Recent experiments show that the conductivity 
in graphene is modified in samples suspended above the substrate [25 1126] . The presence 
of water molecules can screen charged impurities [27] . as well as change the adhesion 
between the graphene layer and the substrate [28ll23] . A detailed investigation of the 
effects of different substrates will be very helpful for the understanding of the carrier 
transport in graphene. 



3 Scattering processes in graphene 

3.1 General framework 

In the following, we analyze different scattering mechanisms which may be present in 
graphene. We first discuss the generalization of the standard partial wave analysis of 
scattering off local potentials in quantum mechanics to Dirac quasiparticles[2^. The 
formalism is rather general, and has already been formulated, using different notation, 
in [30]. Related results can be found in [13|[T4] . We then apply it to different scattering 
processes. The continuum model used neglects scattering between the two inequivalent 
valleys in graphene. We finally give a scheme which takes into account explicitly the 
lattice structure, and which can be used to study scattering processes where intervalley 
scattering is significant. 
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We analyze the phaseshifts induced by a circular potential well in graphene. Using 
cylindrical coordinates, the Hamiltonian in the clean system can be written as: 



ie-'^dr + 

ie'^dr ~ 







'-dj, \ 











-ie-"f'dr - / 

(2) 

where the two first entries correspond to the K point, and the two last ones to the K' 
point, and we are setting the Fermi velocity «p = 1. 

In the following, we study scattering processes which do not induce intervalley 
transitions. Hence, we need only consider one valley. We write the incoming wave as: 



1 



zkr . 



i(n+l)0 -inOk 



(3) 



where the angle 6^. defines the direction of the vector k. In the following, we will set 
6)^=Q. The incoming current is j™ = = 0. 

The outgoing wave is: 



lim 'I'out (r) 

|r| — ^oo 



( fid) 

z7r/4 



2 ^kr 

Tvkr 



E„»"/n[Jn(fcr)+iy„(fcr)le™« 

[/„[J„+i(fcr) +iF„+i(fcr)]e'("+l)« 



(4) 



The outgoing current is jT* = {>I'outWx\<I'out} = 2|/(6»)|2 cos(6'), = {<I'outWy\1'out} 
2|/(6')|^sin(e'). The differential cross section is da{e) = {4\f{e)f )/{nk)de. 
In order to obtain eq.@, we have made the ansatz: 



(5) 



and we have used the expansion: 
lim Jn (kr) 




COS ( kr 



kr 



lim Yn{kr) 



\ / sm I kr 

V Trfcr V 



riTi TV 

~^ 4 



(6) 



We assume that the scattering defect has a finite radius, rg. We can write the 
wavefunctions of an electron with energy vpk outside this radius as a superposition of 
terms: 

[Jn{kr)+RnYn{kr)]e''''' 
i[Jn+i+RuYn+i{kr)] eW)« 



1'nir) 



We now write: 



(7) 
(8) 
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Fig. 1 (Color online). Dependence of the cross section on scattering angle for a circular 
potential well such that k'ro = 1. Red: krg = 0.1 + k'ro. Long dashed (green): kro = 1 + k'ro. 
Short dashed (blue): krg = 10 + k'rg- 



And, using eq.([3j, eq.(|4]) and eq.([7ll, we obtain: 



i - Rn 

, D ^^^^ 



fn 

m = e-^-/'Y.^^ (9) 
Weak scalar potentials satisfy Rn ~ Ri-n- Then, for 6 = tt, we have f{9) = 0. 

3.2 Examples 
3.2.1 Potential well 

We assume that the potential for r < rg is Vg. The wavevector of a state with energy 
E = vpk is k' — {E + Vo)/vp. The matching conditions at r = rg are: 

J„(ferg) + RnYn{kro) = Tn-Mk'ro) 
Jn+l(kro) + RnYn+iikro) = r„J„+i(fcVo) (10) 

and: 

„ ^ Jn{kro)Jn+i(k'ro) - Jn+l (fcrp) Jn(fcVo) 

Yn{kro)Jn+l{k'ro)-Yn+i{kro)Mk'ro) ^ ' 

The largest values of Rn when rg are when n = —1 and n — 0. We expand the 
Bessel functions, and obtain: 

Trfclog (fcrg) 

The cross section is zero for 6* = tt, in agreement with the suppression of backscattering 
in the absence of intervalley transitions. Examples of the cross section due to a circular 
potential are shown in Fig.[T]. 
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Fig. 2 (Color online). Dependence of the cross section on scattering angle for crack with 
zigzag edges. Red: kro = 0.1. Long dashed (green): krg = 1. Short dashed (blue): kro = 10. 



3.2.2 Circular crack 

The boundary conditions for the wavefunction !f'(r) = [!P'i(r), !f'2(r)] at void with zig- 
zag edges, satisfies !i'i(r) = 0. Hence: 

Jnikro) + RnYnikro) = 

Rn = (13) 

Yn(krQ) 

The largest value of Rn is for n — 0. Unlike the previous case, Rq and R-i are not 
equal. To lowest order in rg, the cross section is: 

a{e) ^ 2 (14) 

8fc [log(fero/2) + jf 

where 7 is the Euler-Mascheroni constant. The scattering probability is isotropic. Ex- 
amples of the angular dependence of the cross section for a circular crack are shown in 
Fig. [2]. The behavior of the cross section in this case is the same as that found for a 
strong scatterer in[31j 



3.2.3 Pentagonal cone. 

When inserted into a graphene sheet, a pentagon induces a disclination, and a point 
with finite curvature. The sheet around a pentagon forms a cone. The cone can be 
attached to a fiat sheet with the inclusion of heptagons. We assume that the boundary 
is sufficiently smooth, so that there is not intervalley scattering. The flat surface is 
defined for r > rg, and the conical region for r < rg. Inside the conical region, the wedge 
induced by the dislocation is equivalent to a vortex at the apex with flux — tt/S. 
The matching conditions are: 



Jnikro) + RnYnikro) = TnJn+aikro). 
Jn+likro) + RnYn+likro) = TnJn+l+aikro) 



(15) 



7 



(Tie) 


14 




12 




10 




8 




6 




4 

2 


1 
\ 



n 



3n 



2n 



Fig. 3 (Color online). Dependence of the cross section on scattering angle for conical inclusion. 
Red: fcro = 0.1. Long dashed (green): kro = 1. Short dashed (blue): kro = 10. 



where a = -<P/{2tt) = -1/6. Then: 

^ ^ Jn(kro)Jn+l+a{kro) - Jn+a(fcro)Jn+l(fcro) 

Y„(fcro)J„+i+a(fcro) - Yn+i{kro)Jn+a{kro) 

The cases n and 1 — n are not equivalent, and the scattering cross sections at angles 9 
and 2-K — 6 axe not equal. Examples are shown in Fig. [3]. 



3.3 Weak scatterers. The Born approximation 

3.3.1 General framework 

We expand the outgoing wave as: 

1'out{r)^9„,{v)+ /d2rGo(r-r',t^)V(r>,;„(r') (17) 



where both Gq{v,u)) and V(r) are 2x2 matrices, and uj is the energy of the particle. 
We also have: 

1 [ 2,^^ZS^T + VY^ a^e--A* / 1 
Limr—,ooGr,(r,ui) = — / a k ^ = — - ~ -a 

(18) 

where r — jrj. We analyze the scattering from an angle defined by the unit vector n 
into an angle defined by n. Hence, the vector r is parallel to n. The labels defining 
the coordinates used for the calculation of the Green's function are sketched in Fig. [4] . 
The incoming wavefunction is defined in eq.Q. At long distances, |r| = r ^ cxj, we 
can expand: 

|r — r' I ^ r — r' • n (19) 
where n — r/r. Defining kuj = co/vp, the Green's function becomes: 

G{r~r',.)^ ^-' \/ (20) 

VY\/2%kur \e 1 J 

where 6 is the angle approximately given by the direction of n 



Fig. 4 (Color online). Notation used for the calculation of the Green's function (see text for 
details). 



3.3.2 Scalar potential 

We consider first a potential which does not distinguish the two sublattices: 

V(r) 



V{v) 

y(r) 



The outgoing wave is given by: 



'Pout{r,e) 



r 



\V{v')e-'^ V{v') 



X e 



are ^ ' Viy ) 



1 + e 
e 1 1 + e 



Finally, the amplitude which defines the scattering by an angle Q — is: 



and the cross section 

^{d - din) oc 



are ^ ^ 



V{v') 



1 + e 



y ^2^,^-^/c.(„-n,„)^'^(^,) 



[l + cos(e-e„)] 



and scattering in the backward direction, 9 — = vr, is suppressed. 
The previous calculation can be easily generalized to a potential: 



V(r) = 



V{v) + AViv) 

V'(r) - AV{r) 



We find: 



f{0 - e„,) cx 



ate ^ ' 



+ a r e ^ ' AV(v ) 



I7(r') 



1 - e 



1 + e 



;(e 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



The scattering is isotropic, independent of S — Si„, when the potential modifies only 
one sublattice, V^(r) = ±Z\1/(r). Scattering in the forward direction is suppressed when 
the potential is antisymmetric in the two sublattices, Viy) = 0. 
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3.3.3 Isotropic elastic strains. 

An elastic distortion induces a gauge potential|32): 

N ^ f f3t {Uxx - Uyy + 2iUxy) \ , . 

^ ' \Pt{UxX-Uyy-2iUxy) ) ^ ' 

where u^j are the components of the strain tensor, t is the nearest neighbor hopping 
term in the tight binding hamiltonian, and /3 = dt/dd, where d is the bond length. 
An isotropic ripple with a height profile h{r) leads to (in radial coordinales): 

\f3t{^) e^"f' ; 

This potential is not isotropic, and the scattering does not depend only on the angle 
between the incoming and outgoing directions. We define 

(29) 

The outgoing wave can be written as: 

xe''^"'"""'(,-L) (30) 

we can write: 

(n - ni„)r' = jn - ni„ | r' cos(6l„_„.„ - </>) (31) 



where: 



sin(6l) - sin(6li„) ,/6' + 6l,„, 

tan(6'„_„,„) = — 7Z~^ = -cot (32) 

cos(6^) — cos(6^i„) V 2 ' 



so that: 

+ e,„ TT 



" 2 2 

we can make the change of variables (f) = (f)' -\- [6 -\- ^,j„)/2 — 7r/2 in the integral over 
in ea. pO|l . and obtain: 

r dr dip -===e ' ^ g(r ) y. 

VpV^TTKi^r 
/ gi(20'+20+e,„+7r) _|_ gi(-20'-0-2e„-7r) 
X I gi(20'+e+ei„+7r) _|_ gj(-20'-2e-26li„-7r) 

(r) 

dr' — ;=J2(fci^|n - n'|r').9(''') x 
«FV27rfc(.jr 



(33) 




(34) 
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where J2{x) is a Bessel function. The scattering amplitude can be written as: 



Dp 



and the cross section is: 



f dr'j2{k^\n-n'\r')g{r') (35) 



a(e, 



dr' J2{kuj\n — n'jr')(;(r') 



{l + cos[3(6( + e,„)]} (36) 



This function reflects the threefold symmetry of the honeycomb lattice. 

Near the Dirac point, \ki^\ 0, the cross section in eq.(|36} grows as \kuj\'^. For a 
ripple of height h and size I, we have: 



I r » Z 



(37) 



The hopping t is given approximately, by f ~ up /a, where a is the interatomic distance. 
The total cross section, obtained by integrating eq. (|36p over angles is given by: 



(38) 



^'"'S^'--')' ka « 1 



— TP — 



k^l > 1 



The total cross section has a maximum for ku. 
h ~ Inm and / ~ lOnm, we obtain amax ~ lO/i ~ lOnm 



r\ amax ~ {Pt/vpfh'^/l. For 



3.3.4 Effective magnetic vortex 

Topological lattice defects, such as disclinations and dislocations induce, in addition 
to long range strains, an effective vortex at their core [33ll34ll35j . which mixes the two 
valleys in the case of a disclination. We consider now the scattering by such a vor- 
tex alone, neglecting the effect of the elastic strains. An isotropic distribution of a 
(fictitious) magnetic field induces a potential: 

g{r)e^'t' 
g[r)e-'^ 



nrA)^[ ) (39) 



where g[r) oc (e/c)vprB{r). Using the same scheme as in the previous case, we find: 



I VF\'2TTkur \ g(r)e ^ g{r)e 



,fcc^e^'=- ('■+"'■') ( g{r)e^^^-'^^ g{r)e 



X e 



dr' ^"'L— Ji(fc^|n-nV')./(^') x 
gi(e-ei„+7r)/2 _^ gi(e-ei„+7r)/2 

-i{e+Si„+-ir)/2 _|_ g-i(e+e,„-7r)/2 



(40) 
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Fig. 5 (Color online). Geometry used to calculate scattering amplitudes in the discrete hon- 
eycomb lattice. See text for details. 



SO that: 

f{e - e,^) oc e»(9-e,„)/2 Vfc^ f dr'Mk^ln - n'\r')g{r') (41) 
«F J 

We assume that: 

where / is a dimensionless number of order unity, $o is the flux quantum, and I is the 
redius of the distorted region, we find that the total cross section behaves as: 

kuj f I {kuj^^ kijl <^ 1 
a{k^) ~ { k^fl'^ k^l ~ 1 (43) 



3.4 Lattice effects. 

3.4-1 Green's function formulation 

We can generahze the previous analysis to the discrete honeycomb lattice. The scat- 
tering is no longer isotropic, and only some incident angles can be studied analytically. 
We fix the direction of the incident wave, and use periodic boundary conditions in the 
perpendicular direction. A particular case, where the incident wave is along one of the 
symmetry axes of the honeycomb lattice is sketched in Fig. [5] . The analysis requires 
the calculation of the transmission and refiection coefficients for the different chan- 
nels in the problem. The scheme can easily be generalized to any local lattice defect, 
as sketched in Fig.[S]. This scheme, using the continuum approximation for the local 
Green's function, has been discussed in[3T]. The scattering process can be solved if the 
amplitudes of the reffected and transmitted waves along each of the transverse chan- 
nels of the lattice can be determined. We label the transverse channels by the angular 
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momentum fcj^ . Using periodic boundary conditions with A'' unit cells in the transverse 
direction, the allowed values of k± are fcj^ = 2nn/{Na),n = 0, • • • , A*' — 1; a is the 
lattice constant. 

We analyze next the scattering by a local impurity, and then extend the calculation 
to more complex lattice defects. 



3-4-2 Scattering by an impurity. 

A local impurity perturbs only one site of the lattice. The matching of the incoming 
and the reflected and transverse waves needs only be defined along a line of sites which 
includes the vacancy. We assume that the incoming wave has A;^^ = 0. There are A'' 
reflection and transmission amplitudes, , Tj.^ . The continuity of the wavefunctions 
allows us to deflne A'^ equations: 

l + Ro = To 

Rk^ = Tk^ fc± / (44) 

The other A'^ equations needed to determine uniquely the values of R^^ , Tfc^ are given 
by the condition that the full wavefunction must be an eigenvector with energy e: We 
use ea. (|44p in order to eliminate the R^.^ , and obtain: 

Go\k^,0,e)xTk^ =vJ2Go{k'±,0,€) k^^O (45) 

where V is the strength of the impurity potential, and: 

Go(fei,0,a;) = V i— — (46) 

^ u- e(k»,k^) 
fell 

is the Green's function of the problem, resolved in transverse momentum, and projected 
on the row where the impurity is located. The energies e{k^^,kj_) correspond to the 
unperturbed hamiltonian. 
The solution of ea. (l45|l is: 

r,,^i+ yGoik,,Oe) 

l-^^Efe'^Go(fc^,0,6) 

VGo{k^,0,e) 
1 - V^Z^fc^ Go(fc^,0, e) 

At low energies, we can write: 

e{k^^,k^)^vp^kj+kl (48) 
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and: 




(49) 



where /I is a high momentum cutoff of the order of the inverse of the lattice constant, 
A ~ a-^. 

In order to obtain the total transmitted current as function of the incoming current 
one must multiply iTl-yl'^ by the outgoing current, jj^^^ and divide it by the ingoing 
current, jo- In the tight binding model which describes each transverse channel, we 
have that j|| — G'^^{kj^,Q,Lo). Hence, we can write for the cross section: 



(T(fc|| 



Go(fco,6)y^Go(fc||,6) 



(50) 



and, in order to obtain the dependence of the cross section on the outgoing angle, 
we must take into account that: 



a{e,e)d9 = cr(fc||)dfc|| = a{kpe)^k'^^ + k]_ cos{d)de 



(51) 



Using these expressions, we find that the scattering by a lattice impurity is isotropic 
at low energies, and valley independent (note that the calculation does not distinguish 
the valley index of the transmitted wave) . For V <^t, where t is the nearest neighbor 
hopping, the total cross section is given, approximately, by: 



kf^+klV^ 



(52) 



and, for V t, we find, neglecting logarithmic corrections: 

1 



fc2 



(53) 



In agreement with the results obtained for a crack in the continuum limit, eq. (|14p . 



4 Conclusions 

We have analyzed scattering processes which will affect the mobility of carriers in 
graphene. We show that localized defects can be classified into at least two types with 
opposite dependence of the cross section on density: i) weak scatterers, where the cross 
section grows as the square root of the density, and strong scatterers, where the cross 
section decreases as the square root of the density. The resulting conductivity can be 
either density independent, or grow linearly with density. 

We have also studied defects with internal structure, such as those induced by 
elastic strains or ripples, where the perturbation couples to the Dirac quasiparticles as 
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a) b) 

Fig. 6 (Color online). Geometry used to calculate scattering amplitudes in the discrete hon- 
eycomb lattice in the presence of bond disorder. Left: perturbation of a single bond. Right: 
Stone- Wales defect. See text for details. 



an effective gauge field. We find that for ripples with a characteristic size I, the cross 
section is highest for a density such that fcp^ ~ I. The scattering at each valley does 
not show a symmetry between 9 and —0, where 9 is the incident angle, as expected 
from general symmetry considerations. 

We have finally shown that the classification of short range potentials into weak 
and strong scatterers with different dependence on carrier density remains unchanged, 
even when intervalley scattering is important. 
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